Improving performance of probabilistic programmable quantum processors 
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We present a systematic analysis how one can improve performance of probabilistic programmable 
quantum processors. We generalize a simple Vidal-Masanes-Cirac processor that realizes U(l) ro- 
tations on a qubit with the phase of the rotation encoded in a state of the program register. We 
show how the probability of success of the probabilistic processor can be enhanced by using the 
processor in loops. In addition we show that the same strategy can be utilized for a probabilistic 
implementation of non-unitary transformations on qubits. In addtion, we show that an arbitrary 
SU(2) transformations of qubits can be encoded in program state of a universal programmable 
probabilistic quantum processor. The probability of success of this processor can be enhanced by 
a systematic correction of errors via conditional loops. Finally, we show that all our results can be 
generalized also for qudits. In particular, we show how to implement SU (N) rotations of qudits via 
programmable quantum processor and how the performance of the processor can be enhanced when 
it is used in loops. 



I. INTRODUCTION 

Classical computers are programmable, that is, the 
task they perform is determined by a set of instructions 
that is sent into the machine along with the data to be 
processed. This is a very desirable feature; we do not 
have to build a different circuit every time we want to 
perform a different procedure. It would be useful to be 
able to develop quantum processors with the same prop- 
erty. 

The development of programmable quantum circuits is 
an area that has attracted attention only recently. The 
basic model for these circuits consists of two parts, a data 
register and a program register. There are two inputs, a 
data state, which is sent into the data register, and on 
which an operation is to be performed, and a program 
state, which is sent into the program register, that speci- 
fies the operation. The first result was due to Nielsen and 
Chuang, who showed that a deterministic universal quan- 
tum processor does not exist [I]. The problem is that a 
new dimension must be added to the program space for 
each unitary operator that one wants to be able to per- 
form on the data. A similar situation holds if one stud- 
ies quantum circuits that implement completely-positive, 
trace-preserving maps rather than just unitary operators 
[2,3]. Some families of maps can be implemented with 
a finite program space, for example, the phase damp- 
ing channel, but others, such as the amplitude damp- 
ing channel, require and infinite program space. If one 
drops the requirement that the processor be determin- 
istic, then universal processors become possible [1,4-6]. 
These processors are probabilistic: they sometimes fail, 
but we know when this happens. 

A number of examples of programmable quantum cir- 
cuits have been proposed. One is a quantum "multime- 
ter" that performs unambiguous state discrimination on 



a set of two states, the set being specified by the program 
[7]. There are also devices that evaluate the expectation 
value of an arbitrary operator, the data representing the 
state in which the expectation value is to be evaluated 
and the program state specifying the operator [8,9]. 

In a probabilistic processor, one measures the output 
program state. If the proper result is obtained, the de- 
sired operation has been performed on the data state, 
and if not, then the output of the data register is dis- 
carded. In this kind of a scenario, one wants the prob- 
ability of successfully performing the operation to be as 
close to one as possible. In fact, what one would like, 
is, given a set of operations that one wishes to perform, 
a procedure for systematically increasing the probability 
of successfully performing these operations. 

In the case of one-parameter unitary groups acting 
qubits this was done Preskill [4] and Vidal, Masanes 
and Cirac (VMC) [5]. Vidal, Masanes and Cirac con- 
sidered the one-parameter group of operations given by 
U(a) = cxp(iaer 2 ), for < a < 2ir, and discussed two 
equivalent methods of making the probability of perform- 
ing U (a) arbitrarily close to one. A circuit consisting of 
a single Controlled-NOT (CNOT) gate, with the control 
qubit as the data and the target qubit as the program, 
can successfully perform U(a) with a probability of 1/2. 
If the procedure fails, however, the data qubit, which was 
initially in the state \ip), is left in the state U(—a)\ip). 
What we can now do, is to send this qubit back into the 
same circuit, but with the program state that encodes 
the operation [/(2a). This also has a probability of 1/2 
of succeeding, and increases the total success probabil- 
ity for the two-step procedure to 3/4. Note that our 
program state has increased to two qubits, one for the 
first step and one for the second. We can continue in 
this way simultaneously increasing the success probabil- 
ity and the size of the program state. It is also possible to 
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design more complicated circuits that perform the entire 
procedure at once, i.e. they have a one-qubit data state, 
an TV-qubit program state, and a success probability of 
1 - (1/2)" [5]. 

Here we would like to extend these ideas in a num- 
ber of different directions. First, we shall show that it 
is possible to boost the probability of sets of nonunitary 
operations. It will then be shown how to increase the 
success probability of operations on qudits. Finally, more 
complicated groups of operations will be considered. 



II. OPERATIONS ON QUBITS 

We shall begin by describing the methods developed 
in [4] and [5] in terms of the formalism presented in [6]. 
There, the input data state is in the Hilbert space 7id, 
the program state in the space TC P , and G is the unitary 
operator, acting on the space 'Hd®'Hp, that describes the 
action of the circuit. This operator can be expressed as 



N 



G = A 3k®\j)p P {k\, 

j,fe=0 



where N is the dimension of H p , Ajk is an operator on 
Tid 1 and {\j)\j = 1, . . . N} is an orthonormal basis for the 
program space. The operators Ajk satisfy [6] 
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(2.2) 



where Id is the identity operator on Tid- If the circuit 
acts on the input state \tp) d ® |S) P , we find that 



N 



G(\iP) d ® |H) P ) =^A J -(S)|V) d ® \j) P , 



where 



N 



(2.3) 



(2.4) 



fe=i 



Let us begin by using this formalism, let us look at a 
CNOT gate and the simplest of the circuits discussed in [5] . 
Both the data and program space are two-dimensional, 
and the data space is the control qubit and the program 
space is the target qubit. Expressing the operator for the 
CNOT gate in the form given in Eq. (2.1), and choosing 
the basis {|0), |1)} for the program space, we find that 



A» = |0><0|; A 01 = \l)(l\; 
A 10 = |l)(l|; A n = \0)(0\. 



(2.5) 



We want to use this circuit to perform the operation U (a) 
and this can be done with the program state 



i(a)) = - 7 =(e-|0)+e--|l)). 



(2.6) 



This gives us the output state 

i 

® |S(a)) p ) =Y,M<*M* ® \j) P (2.7) 

3=0 



where the program operators are 



A (a) 



V2 



V2 



ia p-ia 



|0)(0| 
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U(a) 

U(-a). (2.8) 



Therefore, if we measure the output of the program regis- 
ter in the computational basis and obtain |0), then U(a) 
has been carried out on the data state. This occurs with 
a probability of 1/2. 

If we obtain |1) instead of |0) when we measure the 
program register output, then the operation U(—a) has 
been performed on the data state. We can try to correct 
this by sending the state U(— a)\ip)d back into the same 
(2.1) circuit, but with the program state |S(2a)) p . If we mea- 
sure the program output and obtain |0), then the output 
of the data register is 



U{2a)U{-a)\i>) d = U{a)\i,) d , 



(2.9) 



and this happens with a probability of 1/2. This will 
correct the previous error. 

A circuit that does this all at once can be constructed 
from three qubits and two quantum gates [5]. Qubit 1 
is the data qubit, and qubits 2 and 3 are the program 
qubits. The first gate is a CNOT gate with qubit 1 as the 
control and qubit 2 as the target. The second gate is a 
Toffoli gate with qubits 1 and 2 as controls and qubit 3 
as the target. A Toffoli gate does nothing to the control 
bits, and does nothing to the target bit unless both con- 
trol bits are 1, in which case it flips the target bit. If we 
denote the orthonormal program space basis by 



|0) P 
|1>„ 



|0) 2 |0) 3 ; 
|0) 2 |l> 3 ; 



|2> P = |i) 2 |0) 3 ; 
|3) p = |l) 2 |l> 3 , 



(2.10) 



then this circuit can be described by the operators 



A m = |0)(0|; 
A w = 0; 
A 20 = 0; 



A i 
An 
A 2 i 



0; 

|0)(0|; 
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|1)(1|; A 31 =0; 



A 02 
A 12 
A 22 
A 32 



0; 

|0)(0|; 
0; 



Ai 3 

A 23 
A 33 



The program state is now 



,*(3-2j)q 



0; 

|i)(i|; 

0; 

|0)(0|. 
(2.11) 



(2.12) 



1=0 
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At the output of the processor the program register is 
measured in the computational basis, and only if both 
qubits are found to be in the state |1) does the proce- 
dure fail. The overall probability of succeeding is 3/4. 

Now let us go back to the CNOT gate with a single qubit 
program and consider a more general program state 



1 5) =c|0)+ci|l), 
the operators A (E) and Ai(S) are 



A (E) =c |0)(0|+ci|l)(l| 
A 1 (S)=c 1 |0)(0|+c |l)(l| 



(2.13) 



(2-14) 



These operators are not unitary, but they do have the 
property that A (E)Ai(E) = Ai(E)A (E) = cqcJ. The 
output state of this circuit is given by Eq. (2.7), so 
that it can be used to realize, probabilistically, either 
of the nonunitary operators, A (E) or Ai(E). It also 
suggests that we should be able to apply something like 
the Preskill-Vidal-Masanes-Cirac scheme. In particular, 
suppose we are trying to perform the operation 



B(*) = |0)<0|+*|1><1|. 



(2.15) 



If ci = zc , then A (E) is proportional to B{z). We 
send the data state into the processor and then measure 
the program state in the {|0), |1)} basis. If we get we 
have succeeded, but if we get 1 we have instead applied 
A\{E) to the state. If we fail, however, we can try again. 
We now take the output from our first attempt, which is 
Ax(E)\tp)d, and send it into the processor again, but this 
time with the program state 



|S') = 



1 



1 + \z\ 



1/2 



Z 2 |l)). 



(2.16) 



We again measure the program state, and if we find 0, the 
output of the data register is the desired state, Ao(E.)\ip)d- 
If we failed, that is we found 1, we can try yet again, but 
we need to modify the program state every time we re- 
peat the process. 

Rather than performing this procedure sequentially, 
i.e. sending in the input state, seeing if we succeed, and if 
not trying the procedure again with a modified program 
state, we can again do everything at once by enlarging 
the size of the program space. We shall use a slightly dif- 
ferent processor than the one used by Vidal and Cirac. 
It has the same 4-dimensional program space, but the 
operators Ajk are now given by 

A o = |0)(0|; A 01 = |1)(1|; A 02 = 0; A 03 = 0; 

A w = 0; A n = |0>(0|; A 12 = |1)(1|; A 13 = 0; 

A 2Q = 0; A 21 =0- A 22 = |0> <0|; A 23 = |1)<1|; 

A 30 = |l)(l|; A 31 =0; A 32 = 0- A 33 = |0)(0|. 

(2.17) 

The program state is now 



|^)p — Ck\k) p , 



(2.18) 



k=0 



where Ck+i = zck for k = 0, 1, 2, and normalization then 
requires that 



l c o| 



i-N 2 



(2.19) 



The operation of the processor is given by 



G(W d ® |H) P ) =J2A j (Z)W) d ® \j) p , (2.20) 



where 



Aj^) — ^ CkAjk, 



(2.21) 



fe=0 



and the operators Ajk are given in Eq. (2.17). This pro- 
cessor will perform the operation B(z) with a reasonably 
high probability. In order to sec this, we first note that 
Aj(E) = z j A (E) for j = 0, 1, 2,. This implies that 

2 

G(|V) d ® |S) P ) = Ao(SMd ® (E 

3=0 

+A 3 (E)\tP) d (g> |3) p , (2.22) 

and Aq(E) = CqB(z). At the output of the processor we 
measure the program state in the {\j)\j = 0, ... 3} basis, 
and if we get 0,1 or 2, we have carried out the desired 
operation. If \ip) d — a\0) + f3\l), then the probability of 
success depends on the input state and is given by 



P _ i 1 l z l 6 \ i i ., - . z l\- . 



1- z[ 



(\a\ 2 + \z\'\pn 



(2.23) 



If we average this probability over all input states we find 
that 



1 / 1 



2 V 1 - \z\ 



(i + M 2 )- 



(2.24) 
= 1/2, 



As an example, we can consider the case \z\ 2 
which gives us P suc = 0.7. 

This can easily be generalized to an iV-dimensional 
program. The operators Ajk are now given by 



i4 jfc = <y jifc |o)<o| + <y j+ i, fc |i)<i|, 



(2.25) 



where the addition in the second Kronecker delta is done 
modulo N. These operators satisfy Eq. (2.2), so that they 
define a unitary operator. The program state is now 



JV-l 

I 5 ) = c oE zJ I-?)p' 

3=0 



(2.26) 
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where 



1 



co 



l-\z\ 2N ' 
This yields the following output state 



(2.27) 



N-2 



<8> |H) P ) - coB^Md ® ^ ^I^p 

3=0 

+A N _ 1 (Z)\*p d )®\N-l) p , (2.28) 



where 



^ J v-i(S)=co(^- 1 |0)(0| + |l)(l|). 



(2.29) 



The probability of successfully performing B{z) on 
is given by 



P suc = l-||A JV _i(S)V|| 2 



(l-|zn(|qp|z|^) + |/3| 2 ) 
1 I2JV _ ^ 



(2.30) 



When \z\ = 1, this is equal to 1 — (1/N). An examination 
of Psuc shows that it is an increasing function of N. In 
the case that \z\ = 1 it approaches 1 as N — > oo. This is 
no longer true if \z\ ^ 1; if \z\ < 1, we find that the limit 
is 

P S „ C -1-(1-N| 2 )|/3| 2 = ||B(^)V|| 2 , (2.31) 
and if | z \ > 1, the limit is 

Psuc - 1 - (l - j^j = ±- 2 \\B{z)M\ 2 - (2.32) 

Therefore, only in the case that we are implementing a 
unitary operation can this sequence of processors achieve 
a success probability arbitrarily close to 1. 



III. QUDITS 

We now want to see how these arguments can be gen- 
eralized to higher dimensional systems, and, for the sake 
of simplicity, let us start by examining qutrits. The data 
space is now three-dimensional, and let us take for the 
operators Ajk 



A 00 


= |0)(0|; 


An 


= |i)(i|; 


A 02 


= |2)<2|; 




A w 


= 12X21; 


An 


= |0)(0|; 


Al2 


= |i)(i|; 




A 2 a 


= |i)(i|; 


A 2 i 


= |2)(2|; 


A 22 


= |0)(0|. 


(3.1) 



The general program state is 

|S)=c |0)+ci|l)+c 2 |2) ) 
which gives the program operators 



(3.2) 



4(S) = c |0)(0| 
4(S) = c |2)(2| 
A 2 (S)=c |l)(l| 



ci|l)(l|+c 2 |2)(2|; 
Cl |0)(0|+c 2 |l)(l|; 
hci|2)(2|+c 2 |0)(0|. 



The output state is 



3=0 



® \j) P , 



(3.3) 



(3.4) 



so that if we measure in the program space and get j, 
the output state of the data register is Aj(E)\ip)d- 

Suppose we are trying to apply the operator 4) (2) to 
the input data state. The probability of succeeding is 
(^>|4d(")Ao(S)|^). If we fail, however, we can try again, 
and this will increase the total probability of success. To 
see how this works, let us consider an example. Suppose 
that we measured the program register and got 1 instead 
of 0. That means we now have the state 4(S)|V>)d- We 
take this state and put it through the processor again, 
but with a modified program state 



|S')=c[,|0)+ C ' 1 |l)+ C ' 2 |2). 



(3.5) 



Suppose we now measure the output in the program 
space and get 0. If A (E')Ai(E) oc A (S), then we have 
succeeded on our second try. Noting that 



4(S')=c[,|0)(0|+ C ' 1 |l)(l|+ C ' 2 |2)(2|, 



we see that this condition is satisfied if 



aci 
c 2 



Co = 



ac 2 
co 



(3.6) 



(3.7) 



The constant a is chosen so that |S') is normalized. 

What we can conclude from this is that we can, by trial 
and correction, boost the probabilities of implementing 
operators that are diagonal in the basis {|0), |1), |2)}. In 
the case that the operator we are trying to implement is 
unitary, i.e. \cj\ — l/-\/3, then our probability of success 
at each trial is 1/3, so that our probability of success 
after N trials is 1 — (2/3) N . This probability goes to 1 
as N goes to infinity. These conclusions generalize in a 
straightforward way to qudits. 

We now want to explore increasing the probability of 
successfully performing an operation on qudits by in- 
creasing the size of the program space. The data space is 
now of dimension D, and the orthonormal basis spanning 
it is {|0)(j, . . . \D — l)d}- We shall consider a particular 
kind of operation, one that changes the amplitude of one 
of the basis states, and leaves the rest alone (up to overall 
normalization). Suppose the state whose amplitude we 
want to change is \0)d- The operator we want to imple- 
ment is 



B (z) = z\0) pp (0\+X, 



(3.8) 



where 
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D-l 

X=J2\k) PP (k\. (3.9) 

fe=i 

For our processor, we shall choose the operators Aj k , 
where j and k run from to D — 1 to be 

A jk = S jk X + 5 kJ+1 \0) pp (0\, (3.10) 

were all additions are modulo D. The program state 

N-l 

|S) p = co5> fe |fc} P , ( 3 - U ) 

k=0 

where |co| 2 is given by Eq. (2.27), gives us, for < j < 
N-2 

A 3 (E) = c zW (z). (3.12) 

The probability of successfully performing B (z) on the 
data state \tp)d, Psuc, is 

I |2(JV-1) _ 1 
Psuc= ' \ zV N_ 1 ll^ll 2 , (3-13) 

when z\ ^ 1, and it is (N — 1)/^V when \z\ = 1. In the 
limit that N goes to infinity, P suc goes to one if \z\ = 1. 
If \z\ > 1 we have that 

Psuc^ ^\\B i,\\\ (3-14) 
\z\ 

and if \z\ < 1, then 

Psuc^\\BoM\ 2 - (3-15) 

As before, we see that it is only in the case that the 
operation is unitary that the probability goes to one. 

If we want to modify more than one basis vector am- 
plitude, we can apply these processors successively, each 
designed to modify a single amplitude. In the case that 
all of the operations are unitary, this is a Z?-dimcnsional, 
programmable phase gate, whose probability of succeed- 
ing can be made arbitrarily close to one. 

IV. REALIZATION OF SU(2) ROTATIONS 

In the Vidal-Masanes-Cirac model the angle of the 
U(l) rotation that is supposed to be performed on a qubit 
is encoded in a quantum state of the program. The rota- 
tion itself is then applied on the data qubit via the CNOT 
gate that plays the role of a programmable processor. As 
we have discussed above the probability of success of the 
rotation can be enhanced, providing the data qubit is 
processed conditionally in loops. The dynamics of each 
"run" of the processor is conditioned by the result of the 
measurement performed on the program register. 

In what follows we will show that an analogous strat- 
egy can be applied in the case of the SU(2) rotations 



of a qubit, when the parameters (angles) of the SU(2) 
rotations are encoded in the state of the program. In 
our earlier work [6] we have shown an arbitrary singlc- 
qubit unitary transformation can be implemented with 
the probability p = 1/4 by using a quantum informa- 
tion distributor machine (QID) as the processor. The 
QID is a quantum processor with a single data qubit 
and two program qubits. The quantum information dis- 
tribution is realized via a sequence of four CNOT gates, 
such that firstly the data qubit controls the NOT oper- 
ation on the first and the second program qubits and 
then the first and the second program qubits act as 
the control with the data qubit as the target. At the 
end of this process a projective measurement on the 
two program qubits is performed. The measurement 
is performed in the basis: {|0)|+); |0)|->; |1)|+); |1)|-)} 
(where |±) = (|0) ± |1))/V2). The realization of the 
desired transformation is associated with the projection 
onto the vector |0)|+). In what follows we will explicitly 
show how to correct the cases of wrong results, i.e. of 
projections onto one of the vectors |0)|— ), |1)|+), |1)|— ). 

The action of the QID processor is given by relation 
[10,6] 

3 

G = 2^®|H J -)<H J -|, (4.1) 

1=0 

where aj are standard cr-matrices with a^ = I. The basis 
program vectors \Ej) form the standard Bell basis, i.e. 

|So) = -^(|00) + |ll)); |S X > = ^(|01> + |10»; 

|S Z ) = -^(|00)-|11)); |H W ) = ^(|01)-|10)). 

The general program state |E(/J)) P encoding the unitary 
transformation Up = exp(ij2.a) = cos [il + ism^^.a 
(ji = \p\) is given by the expression 

Sill Li 

l s (A% = cos/z|E ) + i (Mx|S x ) + n v \Zy) + H z \^z)) ■ (4.2) 

A* 

Performing the previously mentioned measurement in the 
program basis 1 0+) , 1 0— ) , 1 1+) , 1 1 — ) we obtain the follow- 
ing unitary transformations 



|0> <5 


5|+> 




-UpW) d ; 


|0> <§ 


»R 




■* VzUpa z \ij)) d 


|1>€ 


5|+> 


\i>)d- 


-» a x Upa x \ip) d 


|1>« 




\1>U- 


-» a y Upa y \il))d 



where 

isin^i 

Up = cos^iI-\ (fJ- x a x + fi y a y + [i z a z ) . (4.3) 

M 

To obtain this simple expression we have used the iden- 
tity aja k aj = —a k if k ^ j. All observed outcomes occur 
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with the same probability, p — 1/4. Using the above no- 
tation the action of the QID can be expressed in the form 

|S(/2)) P -> i |j) p j (4.4) 

where vectors {| j) p } form the basis of Tt p associated with 
the realized measurement. The explicit form of the vec- 
tors is presented in following Section where we discuss a 
general solution of SU(N) rotations of qudits. 

We see that each outcome of the measurement indi- 
cates a different unitary transformation has been applied 
to the data. Once we have obtained a specific result we 
can use the same processor again to correct an incorrectly 
transformed data register and consequently improve the 
success probability. In particular, in the case of the result 
j, the new program register needs to encode the correct- 
ing transformation = Upo-jUto-j. The probability of 
implementing the unitary transformation using one con- 
ditioned loop is given as p(l) = \ + 3^ = j^. Us- 
ing more and more conditioned loops the success prob- 
ability is given by p(n) = YT ]= i 37 = l^jij) 3 = 
3 1 ~i% 4 ^ = ^ — (3/4)™ converges to unity, i.e. p(n) — > 1 
as the number of conditioned loops n goes to infinity. For 
instance, thirty conditioned loops result in the negligible 
probability of failure, p ~ 10~ 4 . 

The example of Vidal, Masanes and Cirac shows us 
that we are able to replace the feedback scenario with a 
probabilistic scenario by using different processors. An 
open problem is whether the same replacement can be 
done in general, or at least for the case of the QID. 

V. SU (N) ROTATIONS OF QUDITS 

In what follow we will show that one can utilize the QID 
for a probabilistic implementation of SU(N) rotations of 
qudits. We start our discussion with a brief description 
of the QID in the case of qudits. First, we introduce a 
generalization of the two-qubit CNOT gate [10] for qudits. 



This is a conditional shift operator defined with a control 
qudit "a" and the target qudit "6" 

N-l 

D ab = \k)a(k\®\(m + k)modN) b (m\ , (5.1) 

k,m— 

which implies that 

N-l 

D L= \k)a(k\®\(m- k)mod N) b (m\. (5.2) 

k,m— 

^From this definition it follows that the operator D ab acts 
on the basis vectors of a qudit as 

D ab \k)\m) = \k)\(k + m)modN) , (5.3) 

which means that this operator has the same action as 
the conditional adder and can be performed with the help 
of the simple quantum network discussed in [11]. Note 
that for N > 2 the two operators D and differ; they 
describe conditional shifts in opposite directions. There- 
fore the generalizations of the CNOT operator to higher 
dimensions are just conditional shifts. 

Following our earlier work [6,10] we can assume the 
network for the probabilistic universal quantum proces- 
sor to be 

DnD^DmDu . (5.4) 

The data register consists of system 1 and the program 
register of systems 2 and 3. The state |Sy)23 acts as the 
"software" that caries the information about the opera- 
tion V to be implemented on the qudit data state \^>}i- 
The output state of the three qudit system, after the four 
controlled shifts are applied, reads 

|fi>i23 - £> 3 i J D| 1 £» 13 £>i2|*)i|S y ) 2 3 . (5.5) 

The sequence of four operators acting on the basis vec- 
tors gives |n)i|m) 2 |fc) 3 as 



D 31 D\ l D 13 D 12 \n) 1 \m) 2 \k) 3 = |(n - m + k)mod N)i \(m + n)modN) 2 \(k + n)modN} 3 . 



(5.6) 



r 



We now turn to the fundamental program states. A ba- 
sis consisting of maximally entangled two-particle states 
(the analogue of the Bell basis for spin-i particles) is 
given by 

1 N ~ 1 2ir 
\E mn ) = — exp(i-^mfc)|fc)|(fc-n)modAr), (5.7) 

^ fc=0 

where m, n = 0, . . . , N — 1. If |S mn ) p is the initial state 
of the program register, and |\f') = Ylj a i\j)<l (here-; as 



usual, J2j \ a :\ 2 = 1) i s tne initial state of the data regis- 
ter, then follows that 



fl23|*)l|-mn)23 

Ea, \2nikm 
v! exp 



jk 



N 



Pi2s\j)\k)\k-n) 



Eaij 2-ikm 
^=exp — 1 ^—\3 - n)\k + 3)\k + 3 - n) 



jk 



N 
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E—2mjm, . 
a 3 ex P Jj \J - »>Fmn> 



JV-l 



jk 



23- 



(5.13) 



= (f/(™>«)|*))|S m „), 
where we have introduced the notation 



m,n— 



(5.8 



JV-l 



f /(™-)=^exp^=| S -n)( S |. (5.9) 

This result is similar to the one we found in the case of a 
single qubit (see previous section). The operators {/( m ' n ) 
satisfy the orthogonality relation 



Application of the processor to the input state |^')i|^y)23 
yields the output state 

I^)i23 = £dm„tf (m ' n) |*>i ® |S mn ) 23 . (5.14) 



Now let us perform a measurement of the program out- 
put in the basis 



Tr 



fjj{m' ,n')^\jj{m,n 



(5.10) 



The space of linear operators T(H) defined on some 
Hilbert space Tt with the scalar product given by (5.10) 
we know as Hilbert- Schmidt space. Thus the unitary op- 



1 N ~ 1 

\®rs)=K E 6XP 



m,n— 



2iri 



(mr — ns) 



N 



|S mn >. (5.15) 



The orthogonality of this measurement basis directly fol- 



erators ^™>")form an orthogonal basis in it and any lowg from ^ orthogonality of the entang i ed basis | Smn) . 



operator V G T(7{) can be expressed in terms of them 



JV-l 



We should also note, that the vectors |$ rs ) itself can be 
rewritten in a factorized form, i.e. 



V= d mn U^ n \ 



(5.11) 



m,n— 



The orthogonality relation allows us to find the expan- 
sion coefficients in terms of the operators 



\§rs) 



r) ® ^ exp 27ri-^ \n-r), (5.16) 

V ^ n=0 



N 



dn 



N' 



-Tr 



V 



Therefore, the program vector that implements the op- 
erator V is given by 



which means that the measurement can be performed 
(5.12) independently on two program qudits. 

In order to clarify the role of the measurement we will 
rewrite the output state of the QID using the basis |$ rs ) 
for program qudits: 

J 



JV-l 

12 rfm,nf/ (m ' n) |*)l|: 


=, mn)23 




m.n—0 






N-l 


N E ex P 

r,s— 




E <wf/ (m,n) i*>i 

m.n—0 


-2m 



. (mr — ns) 



I* 



rs/23 



JV-l JV-l 



^ EE { cx p 



Taking into account that 



r,s — m,n— 



,i t 



. (mr — ns) 



N 



N 

rf m ,nf/ (m ' n) }|*)l|^ S )23 



jjfog) [/(m,«)/j(P-9) = exp 

and choosing p = s and q = r we find 



. (mg — np) 



cxp 



iV 



-2ttz 



(mr — ns) 
iV 



Finally, the output of the QID can be rewritten in the form 

AT-l r 



- IV — J. r- ,■ 

Pl2 3 |*)l|Sy) 23 = ^ E ^ ( '' r) ^ ( C/(5 ' r) ) 



r.s— 



|^>l|*r- S >23 , 



(5.17) 



(5.18) 



(5.19) 



(5.20) 
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from which it is clear that if the result of the measure- 
ment of the two program qudits is | < E> rs ) 2 3, then the sys- 
tem (data) is left in the state \lj( s ^V (U (s ^} 

Obviously, if s = r = 0, then the operator V is applied 
on the data qudit. The probability of this outcome is 
1/N 2 . For all other results of the measurement the data 
qudit is left in the state given above. One can use these 
output states with a modified program state to improve 
the performance of the programmable processor. Specif- 
ically we have to use the new program state |Sy' s ^} that 
is chosen after taking into account the result of the pre- 
vious measurement. This program state has first to "cor- 
rect" the wrong realization of the operation V during the 
previous "run" of the processor and then apply (proba- 
bilistically), the original operation V. For this reason, 
the new program state has to perform the operation 

-1 



y(r,s) _ y 



(5.21) 



This process of error correction (conditional loops) can 
be used K times and the technique of conditioned loops 
can be exploited in order to amplify the probability of 
success. Applying the processor K times the probability 
of a successful application of the desired SU(N) operation 
V reads p(K) = 1 - (1 - l/N 2 ) K . 



of qubits using the probabilistic programmable processor 
with the quantum program register initially prepared in 
states that carry the information about the operation 
to be performed on the data. It is also possible to use 
enlarged programs to increase the probability of success 
without the use of loops. In this case the measurement 
performed on the program register has to be modified 
accordingly. We have shown that if the processor is used 
in loops with properly chosen program states one can 
improve the performance of the quantum programmable 
processor so that the probability of failure decreases ex- 
ponentially with the number of program qudits that store 
the information about transformation on the data qudit. 
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VI. CONCLUSIONS 

In this paper we have analyzed a probabilistic pro- 
grammable quantum processor. We have shown how to 
encode information about the quantum dynamics V to 
be performed on a quantum system (data register) in 
the state of another quantum system (program register). 
This information is stored in such a way that the program 
can be used to probabilistically perform the stored trans- 
formation on the data. In our paper we have analyzed 
systematically how to perform U(l) rotations of qubits 
and qudits and one-parameter families of nonunitary op- 
erations when the angle of rotation is encoded in states of 
quantum programs. In addition we have shown how to in- 
crease the probability of success when the quantum pro- 
cessor is used in loops with updated program states. We 
have generalized the whole problem and we have shown 
that one can use a very simple quantum processor, the so 
called quantum information distributor, to perform arbi- 
trary SU(2) rotations of qubits as well as SU(N) rotations 
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